Kolmogorov and Bolgiano scaling in thermal convection: 
the case of Ray leigh- Taylor turbulence 



O 



m 

Q 
U 



> 
in 



G. Boffetta\ F. De Lillo\ A. Mazzino^, S. Musacchio^ 
^ Dipartimento di Fisica Generale and INFN, Universitd di Torino, via P.Giuria 1, 10125 Torino (Italy) 
^ Dipartimento di Fisica, Universitd di Genova, INFN and CNISM, via Dodecaneso 33, 16146 Genova (Italy) 
'•^^ CNRS, Lab. J. A. Dieudonne UMR 6621, Pare Valrose, 06108 Nice (France) 

We investigate tiie statistical properties of Rayleigh- Taylor turbulence in a convective cell of high 
aspect ratio, in which one transverse side is much smaller that the others. We show that the scale 
of confinement determines the Bolgiano scale of the system, which in the late stage of the evolution 
is characterized by the Kolmogorov-Obukhov and the Bolgiano-Obukhov phenomenology at small 
and large scales, respectively. The coexistence of these regimes is associated to a three to two- 
dimensional transition of the system which occurs when the width of the turbulent mixing layer 
becomes larger that the scale of confinement. 



Turbulent thermal convection appears in many natu- 
ral phenomena, from heat transport in stars to turbu- 
lent mixing in the atmosphere and the oceans, and in 
technological applications [IHll- Turbulent convection is 
driven by buoyancy forces generated by temperature fluc- 
tuations. These are then mixed by the turbulent flow 
itself up to small scales at which molecular diffusivity 
becomes important. A fundamental problem in thermal 
convection is the determination of the statistical proper- 
ties of velocity and temperature fluctuations in the iner- 
tial range of scales in which turbulent mixing is at work. 

A first step in this direction was done by Obukhov 
[3| and Corrsin [5] who generalized the Kolmogorov ar- 
gument for the statistics of a temperature field in the 
so-called passive limit, in which the effects of the buoy- 
ancy forces on the velocity field are neglected [gK An 
alternativeprediction was proposed by Bolgiano [7] and 
Obukhov '&\ , in discussing the statistics of velocity and 
temperature fluctuations in a stably stratifled atmo- 
sphere. The buoyancy forces allow to introduce in the 
inertial range a characteristic scale, the Bolgiano scale 
Lb, above which the statistics of the velocity and temper- 
ature is determined by the balance between the buoyancy 
and inertia forces. In spite of many years of experimen- 
tal and numerical investigations, no clear evidence of this 
scenario has been given [9|. 

In this Letter we show that in three-dimensional 
Rayleigh- Taylor turbulent convection the Bolgiano scale 
is determined by the geometrical scales of the setup. By 
confining the flow in a box with one dimension (e.g. y) 
much smaller than the other two, the scale Ly becomes 
the Bolgiano scale of the system. By means of state of the 
art numerical simulations of primitive equations we find 
coexistent Kolmogorov-Obukhov scaling at scales smaller 




than Ly and Bolgiano scaling at scales larger than Ly. 
Our geometrical interpretation of the Bolgiano scale sug- 
gests a new direction for numerical and experimental in- 
vestigations of scaling properties in thermal convection. 

Rayleigh- Taylor (RT) turbulence is one of the simplest 
configurations of thermal convection in which a cold. 



FIG. 1: Central part of a vertical {x,z) section of the tem- 
perature field (white=hot, black=cold) from a direct numer- 
ical simulations of Boussinesq equations ^ at resolution 
4096 X 128 X 8192 taken at time t = S.Or (see text for the def- 
inition of r). Parameters are Pg = 0.5, 1/ = K = 2 X 10"^ and 
the initial temperature jump is = 1. The picture clearly 
shows that large-scale 2D structures coexist with small-scale 
3D turbulence. The segment represents the transverse size 

Ly . 



heavier layer of fluid is placed on the top of an hot, lighter 
layer in a gravitational field. Rayleigh- Taylor instability 
occurs in several phenomena ranging from geop hysics, to 
astrophysics to technological applications |10l4l2| . The 
gravitational instability develops in an intermediate layer 
of turbulent fluid (the mixing layer) the width of which 
grows in time. 

We consider miscible RT turbulence at low Atwood 
numbers. Within the Boussinesq approximation, the 
equations for the dynamics of the velocity fleld u coupled 
to the temperature field T(x, t) (which is proportional to 
the density p via the thermal expansion coefficient /3 as 
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FIG. 2: Temporal evolution of the mixing layer width h{t) 
(blue, dashed line) and the ratio of y component to x compo- 
nent of rms velocity fluctuations (red, solid line). Inset: Ratio 
between kinetic energy growth rate dE/dt and energy dissipa- 
tion rate e as a function of time. The horizontal line indicates 
a constant ratio ~ 0.75 and the continuous line represents the 
scaling t^^^ . 



FIG. 3: Third order longitudinal velocity structure function 
Ss,{r) (red circles) and —S^ir) (blue squares) computed in 
the mixing layer. The two lines represent Kolmogorov scaling 
r (blue dotted) and Bolgiano scaling r^^^ (red continuous). 
Inset: contribution to the energy flux in Fourier space n(fc) 
from the nonlinear term (red, continuous line) and from the 
buoyancy term (blue, dashed line). 



p = po[l ~ P{T — To)] where po and Tq are reference 
values) read: 

dtVL + xiVu = -(l/po)Vp + i/V2u-/3gT 

dtT + M-VT = kV^T (1) 

together with the incompressibility condition V • u = 0. 
In ([1]) g = (0, 0, —g) is gravity acceleration, v is the kine- 
matic viscosity and k is the thermal difFusivity. The ini- 
tial condition for the RT problem is an unstable temper- 
ature jump r(x. 0) = —{6{)/2)sgn{z) in a fluid at rest 
u(x,0) = 0. 

As the system evolves, the available potential energy 
P = —/3g{zT) is converted into kinetic energy E = 
(1/2) (|up) at a rate that can be estimated from the en- 
ergy balance: 



dP . , ^> dE 
— = Pg{wT) = — 



(2) 



where w is the vertical velocity and e = t^HdaUp)'^) is 
the viscous energy dissipation rate. From the dimen- 
sional balance between the loss of potential energy and 
the increase of kinetic energy one has that typical veloc- 
ity fluctuations grow as Urms — Pg^oti and therefore the 
width of the turbulent mixing layer h{t), shown in Fig. [21 
grows following the accelerated law h{t) ~ fSgOoP. The 
integral scale L{t) of the turbulent flow, defined as the 
largest scale on which kinetic energy is injected, is ex- 
pected to grow proportionally to the geometrical scale 
h{t) [13. 

According to the phenomenological theory of RT tur- 
bulence, developed in 14| , the scaling behavior of the ve- 
locity and temperature fluctuations in the range of scales 



between the integral scale L and the dissipative scale ij 
strongly depends on the dimensionality of the flow. 

For the three-dimensional (3D) case, one assumes that 
the buoyancy force /3gT balances the inertia term in ([1]) 
at the integral scale L{t) and becomes negligible as the 
cascade proceeds towards small scales, consistently with 
the Kolmogorov-Obukhov phenomenology. For velocity 
and temperature fluctuations Su{r) = u{x + r) — u{x) (u 
denoting one velocity component) and 5T{r) — T{x + 
r) — T{x) one therefore expects [l4| 
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6u{r) ~ "r 
ST{r) 



(3) 



where the energy dissipation rates e grows in time as 
e{t) ~ {l3g6o)'^t, following adiabatically the dynamics of 
the large eddies, while the tem per ature dissipation rates 
decreases as exit) ~ 0ot~^ fuui^. 

This scenario is not consistent in two dimensions (2D) 
where kinetic energy is transferred toward large scales 
developing an inverse cascade In this case the 

buoyancy term injects energy at all scales generating a 
non-constant-in-wavenumber energy flux. As a conse- 
quence, velocity and temperature fluctuations follow the 
Bolgiano-Obukhov scaling 14 1 



6u(r) '. 
5T{r) 
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(4) 



which has been verified in numerical simulations of 2D 
RT turbulence 



Let us consider now a convective cell with large aspect 
ratio Lx,Lz ^ Ly. At short times, when Lit) < Ly 
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the turbulent flow in the mixing layer can be considered 
three-dimensional and a direct cascade with Kolmogorov- 
Obukhov scaling ^ is expected. At later times, when 
L{t) > Ly the flow cannot be fully three dimensional 
at large scales, as a consequence of the geometrical con- 
strain in the y direction. The question arises whether 
the two-dimensional phenomenology appears at scales 
Ly < r < L{t). For this to happen most of the injected by 
buoyancy forces should go to large scales producing an in- 
verse cascade with Bolgiano-Obukhov scaling (|4]). How- 
ever, a residual direct cascade to scales r < Ly should 
be present with a flux e{t) given by the matching of the 
scaling ([3]) and (|4]) a.t r = Ly-. 

e{t) ^ iPgeoP'L^y/'r^^' (5) 

The time of the transition from 3D to 2D behavior is 
given by continuity requirement in energy dissipation, 
i.e. equating ([5]) with 3D dissipation [PgO^Yt which gives 
T = {Lyl{aPgeo)Yl'^ where a ~ 0.02 is a dimensionless 
number obtained from numerical simulations [15|]. 

Summarizing, the long-time behavior of RT turbulence 
with large aspect ratio is the following. A small scales 
T] < r < Ly a. three-dimensional direct cascade with 
Kolmogorov-Obukhov scaling is expected. At large scales 
Ly < r < L{t) ^ a two-dimensional inverse cascade 
with Bolgiano-Obukhov scaling should be observed. 

The above predictions have been tested against the 
results of state of the art, high resolution numerical sim- 
ulations of equations ((1]), is a large aspect ratio geometry 
with Lz/Lx = 2, Ly/Lx = 1/32, in which the flow thus 
results strongly confined in the y direction. The integra- 
tion of equations (HI, discretized on a 4096 x 128 x 8192 
grid with periodic boundary conditions, has been per- 
formed with a fully parallel pseudospectral code, with 
2/3-dealiasing, running on a IBM-SP6 supercomputer. 

Figure [T] shows a vertical section of the temperature 
field in the late stage of the simulations. Large scales, 
2D structures are clearly observed. 

As shown in Figure [2l at t ~ r, when the mixing layer 
scale h{t) becomes of the order of Ly the fiow becomes 
increasingly anisotropic with Uy <C Ux,Uz. In this con- 
ditions, we observe a transition from 3D to 2D turbu- 
lent behavior, clearly signaled by a change in the ratio 
between the energy growth rate dE/dt and the viscous 
dissipation rate e (the energy flux to small scales in the 
direct cascade). In the 3D regime both these quantities 
grow linearly in time and therefore their ratio is con- 
stant. After the transition the inverse cascade sets in 
and {dE/dt)/e ~ as follows from ([2]) and ([5]). Both 
behaviors are evident in Fig. [2] 

In the late stage of the evolution, when L{t) > Ly 
we expect the simultaneous presence of a direct and an 
inverse cascade in the two range of scales r < Ly and 
r > Ly, respectively. This can be verified by com- 
puting the scale dependent energy flux, given by the 
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FIG. 4: Second-order longitudinal velocity structure function 
S2{r). The two lines represent Kolmogorov scaling r^^'^ (blue 
dotted) and Bolgiano scaling r"''^ (red continuous). Upper in- 
set: Second-order structure function 52 (r) compensated with 
Kolmogorov (red crosses) and Bolgiano (blue squares) scal- 
ing. Lower inset: Fourth-order structure function 54 (r) com- 
pensated with Kolmogorov (red crosses) and Bolgiano (blue 
squares) scaling. 

third-order structure function of longitudinal velocity in- 
crements (i.e. taken along the local velocity direction) 
S'3(r) = (((5|[u(r))^). For isotropic three dimensional tur- 
bulence, the classical result due to Kolmogorov predicts 
53(r) = -(4/5)£r. 

As shown in Fig. [31 at small scales r < Ly 5*3 (r) is 
negative and, in a narrow range of scales, compatible 
with the Kolmogorov prediction S3{r) r. At scales 
r > Ly, 5*3 (r) becomes positive, signaling the reversal 
of the energy cascade, and displays a scaling behavior 
Ssir) r^/^ consistent with 

The inset of Fig. |3] shows the contributions of the in- 
ertia and buoyancy terms to the energy flux in Fourier 
space n(fc) = (d/dt) E{p)dp where E{p) is the energy 
spectrum and time derivative is computed by taking into 
account, separately, the non-linear and buoyancy terms 
of ([T]). Al low wavenumbers kLy < 1 the buoyancy contri- 
bution is dominant, and the negative sign of the inertial 
contribution to the energy flux confirms the presence of 
a 2D inverse cascade. At high wavenumbers kLy > 1 the 
buoyancy contribution becomes sub-dominant, and one 
recovers a constant positive flux characteristic of the 3D 
regime [isj . 

The coexistence of Kolmogorov-Obukhov scaling at 
small scales and Bolgiano-Obukhov scaling at large scales 
is confirmed by behavior of the structure functions of lon- 
gitudinal velocity increments Sp{r) = (((5||M(r))P) and of 
temperature increments Sp{r) = {{6T{r))P), shown in 
Figures m and [5] The transition between the two regimes 
occurs at the Bolgiano scale which is found to be r ~ Lj,. 

In three dimensional turbulence, small deviations from 
the dimensional predictions are expected in the scaling of 
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FIG. 5: Second-order temperature structure function S2{r) 



The two lines represent Kolmogorov scaling 



r2/3 



(blue dot- 



ted) and Bolgiano scaling r^''^ (red continuous). Lower inset: 
Second-order structure function S2{r) compensated by Kol- 
mogorov (red crosses) and Bolgiano (blue squares) scaling. 
Upper inset; Fourth-order temperature structure function to- 
gether with power laws corresponding to best fit scaling at 



small scales r 
continuous) . 



(blue dotted) and at large scales 



(red 



the velocity fluctuations, indicating the presence of small 
scale intermittency [l8| . Here, we found much stronger 
corrections in the statistics of temperature fluctuations, 
whose 4-th order structure function strongly differs from 
the dimensional scaling at small scales with a best fit 
exponent 0.9 close to the corresponding exponent for 
passive scalar in three-dimensional turbulence. At large 
scale, temperature structure functions show strong fluc- 
tuations due to the presence of regions of unmixed fluid 
within the mixing layer, as shown in Fig.[T] Nonetheless, 
a very short range of scaling also for the 4-th order struc- 
ture function is observed with a scaling exponent close to 
the intermittent value 0.6 < 4/5 measured in pure 2D RT 
simulations [l7l |. 

Our numerical findings supports the phenomenological 
prediction that in RT convection the Bolgiano scale is de- 
termined by the aspect ratio of the convective cell. The 
presence of Bolgiano-Obukhov scaling is associated to a 
dimensional transition of the flow which occurs when the 
width of the mixing layer becomes larger than the confin- 

3 



ing scale Ly 



20j . This poses the intriguing question 
on whether in generic convective systems the Bolgiano- 
Obukhov phenomenology could be observed whenever 
the turbulent flow is confined by geometrical constraints 
and/or physical mechanisms (such as rotation or strati- 
fication) in convective cells with small aspect ratio. In- 
deed, recent experiments in soap film convection observe 
Bol gian o scaling for large values of the Rayleigh number 
2lL |22|. Despite the similarities with our results, we re- 
mark the presence of important differences, as the exper- 
imental data show both the presence of intermittency in 



the velocity field 2l| and the absence of intermittency in 
the temperature field 22|. On the contrary, our simula- 
tions shows the absence of intermittency for the velocity 
field which performs an inverse cascade and are compat- 
ible with some intermittency for the temperature field 
in the direct cascade. It would be therefore extremely 
interesting to have a direct comparison of experimental 
and numerical data on thermal convection in quasi-two- 
dimensional flow which would give new insights on the 
fundamental issue of intermittency in turbulent convec- 
tion. 

Numerical simulations have been done at the Cineca 
under the IIPC-EUR0PA2 project (project number: 
228398) with the support of the European Commission - 
Capacities Area - Research Infrastructures. We thank G. 
Erbacci and the staff at Cineca Supercomputing Center 
(Bologna, Italy) for their support. 
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